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' Abstract 

^ ' We study the finite temperature free energy and fermion number for 

, Dirac fields in a one-dimensional spatial segment, under two different 

members of the family of local boundary conditions denning a self-adjoint 
Euclidean Dirac operator in two dimensions. For one of such bound- 
^ . ary conditions, compatible with the presence of a spectral asymmetry, 

^j- ■ we discuss in detail the contribution of this part of the spectrum to the 

(*~) zeta-regularized determinant of the Dirac operator and, thus, to the finite 

temperature properties of the theory. 
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1 Introduction 



> 

X 

When the Euclidean Dirac operator is considered on even-dimensional com- 
pact manifolds with boundary, its domain can be determined through a family 
of local boundary conditions which define a self-adjoint boundary problem 1 
(the particular case of two-dimensional manifolds was first studied in (2. ) . The 
whole family is characterized by a real parameter 9, which can be interpreted 
as an analytic continuation of the well known 8 parameter in gauge theories. 
These boundary conditions can be considered to be the natural counterpart in 
Euclidean space of the well known chiral bag boundary conditions. 

Recently, it was shown |3j that the boundary problem so denned is not only 
self-adjoint, but also strongly elliptic @J in any even dimension. Also in 
reference [3], the meromorphic properties of the associated zeta function were 
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determined for manifolds of the product type. For particular non-product man- 
ifolds, heat kernel coefficients and zeta functions were treated in [Hj. Anomalies 
were also studied recently in [7]. 

One salient characteristic of these local boundary conditions is the generation 
of an asymmetry in the spectrum of the Dirac operator. For the particular case 
of two-dimensional product manifolds, such asymmetry was shown, in [jJJ, to be 
determined by the asymmetry of the boundary spectrum. 

The aim of this paper is twofold: 

i. To discuss a physical application of Euclidean bag boundary conditions in 
two dimensions, with emphasis on the effect of the spectral asymmetry. 

ii. To (partially) answer the question posed in [5], as to whether the fermion 
number is modified by temperature in low dimensional bags. 

In sectional we present the zero-temperature problem of Dirac fermions sub- 
ject to local boundary conditions with arbitrary values of in 1 + 1-dimensional 
Minkowski space-time, and evaluate the vacuum energy and fermion number 
(equivalently, U(l) charge). This model has recently been considered, in the 
context of brane theory, as the fermionic sector of an N — 2 supersymmetric 
sigma model, coupled to a magnetic field at the ends of the string JU| . 

In section |31 we determine the spectrum of the Euclidean Dirac operator at 
finite temperature for two particular values of 9. With these spectra at hand, 
we perform, in section the calculation of the free energy via zeta function 
regularization. 

Section 03 is devoted to the evaluation of the free energy in the grand- 
canonical ensemble and the fermion number, through the introduction of a 
chemical potential (The finite-temperature fermion number for a different lo- 
cal boundary condition, not leading to a self-adjoint operator, was treated in 

EH)- 

Finally, section [HI presents a discussion of the main results. 



2 Definition of the problem in Minkowski space- 
time 

We will use the metric (-,+), and choose for the Dirac matrices 

lM=i°i, 1m = -P2 and 7^=7^7^ = 0-3. (2.1) 
The action, for Dirac fermions coupled to a background field is given by 

S M = i J d 2 x^{ip- 4)* . (2.2) 

Let us first particularize to the free case. The Hamiltonian can be deter- 
mined from the classical equation of motion iftSf = 0, by proposing ^(x , x 1 ) = 



2 



e tEx ip( xl ): with < x 1 < L. Thus, one gets the Hamiltonian H — ijMdi. 

(2.3) 



Its eigenfunctions are of the form 



Ae~ ltjX 



Be 1 



The boundary conditions will be taken to be 
1 



pijhddo.L 



(2.4) 



Now, it is easy to see that the eigenfunctions of the Hamiltonian depend only 
on the difference 8 = 8l — 8q, since the overall phase can always be eliminated 
through a constant chiral transformation. We will thus consider 



;(l-i72f)tf 







= 0. 



Once such boundary conditions are imposed, the eigenfunctions read 



-A„e iE " xl 



where 



E., 



Tin 

T 



2L 



with 



-co oo . 



(2.5) 



(2.6) 



(2.7) 



To evaluate the vacuum properties at zero temperature, let us first consider 
< 8 < 2ix. Then, when defined through zeta function regularization |12| . the 
Casimir energy is given by |13j 



E c m 





n 




s=-l 






-1,1 


e 

~ ~2~t~t 


-( 


8 2 


9 






2tt 2 


- 4 

TT 


\)- 



2tt 



(2.8) 



where Ch{s, u) is the Hurwitz zeta function |14) 

In particular, Ec(tt) = — 23l ■ Note that, for —2tt < 8 < 0, the replacement 
8 — > 2tt + 8 must be performed. For 8 — 



(2.9) 
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where is the Riemann zeta function. So, the Casimir energy is continuous at 
9 = 0. 

As for the vacuum expectation value of the fermion number (C7 (1) charge), 
again for < 8 < 2tt, one has JS] 



N(9) = 




(2.10) 



Also in this case, for — 2ir < 8 < the replacement 8 — > 27r + must be 
performed. For 8 = it, one has N(ir) = 0. But, at variance with the Casimir 
energy, the fermion number is discontinuous at 8 = 0. In fact, in this case, 
apart from a symmetric nonvanishing spectrum, a zero mode of the Hamiltonian 
appears, which is its own charge conjugate. As a consequence, iV(0) = ±| |T5] . 

In what follows, we will concentrate on two values of 8, i.e., 8 = and 8 = n, 
to study the effect of a nonvanishing temperature on both vacuum quantities. In 
both cases, the Euclidean Dirac operator will turn to be self-adjoint, as shown in 
X . Note these two boundary conditions are the ones corresponding to Ramond 
and Neveu-Schwarz strings |1U) . 



3 Finite temperature. Spectrum of the Dirac 
operator 

In order to study the effect of temperature, we go to Euclidean space, with 
the metric (+,+). To this end, we take the Euclidean gamma matrices to be 
7o = — = <J\, 71 = — 7m = °2- Thus, the Euclidean action is 

S E = -iS M = J d 2 x 4)* , (3.1) 

where the Euclidean zero component of the gauge potential (Aq) is related to 
the corresponding Minkowski one by Aq = —iA^ 1 . 

We start by treating the free case; then, the partition function in the canon- 
ical ensemble is given by (see, for example, |2(J| and references therein) 

log Z = log det(ip BC ). (3.2) 

Here, BC stands for antiperiodic boundary conditions in the "time" direction 
(0 < xq < P, with f3 = ^) and, in the "space" direction (0 < X\ < L). 

^(l+7o)* 



4 



5(!±7o)* 



0. (3.3) 



In the last equation, the plus sign corresponds to the case 9 = 0, while the 
minus sign corresponds to 8 = n. 

In order to evaluate the partition function in the zeta regularization ap- 
proach, we first determine the eigcnfunctions, and the corresponding eigenvalues 
(uj), of the Dirac operator 

m= f id Q + d 1 \f V {x Q ,x 1 )\ =u f^x l )\ 
\ido-di J \ x{%o,%i) ) V Xl^o,^iJ / 

To satisfy antiperiodic boundary conditions in the xq direction, we expand 
*(x ,x 1 ) = J2e iXxo iP(xi), (3-5) 

A 

with 

A/ = (2Z + 1)^, I = -oo,...,oo. (3.6) 

After doing so we have, for each Aj, 

(-Xi+di)x = utp 

(-Xt-9i)<P = wx- (3-7) 

3.1 # = 

It is easy to see that, with this boundary condition, no zero mode appears. For 
uj ^ one has, from |3.7|>. 

9^ = —K 2 (p 

X=--(k+dx)<p, (3.8) 

where k 2 = uj 2 - Xf. 

For k / 0, one has for the eigenvalues 



(nn\ 2 2 

w„,i=±Y^— J + A/ , with n = l,...,oo. (3.9) 

This part of the spectrum is symmetric. In the case k = one has a set of 
Xi-independent eigenfunctions, corresponding to 

oji=Xi. (3.10) 

It is to be noted that u>i = — A; are not eigenvalues; so, this last portion of 
the spectrum will be asymmetric or not, depending on whether the boundary 
spectrum ({A;}) is so. This was to be expected from our result in [H], where 



5 



we proved that, in this case, the contribution of each boundary to the asymme- 
try equals one half the boundary asymmetry, and the contributions from both 
boundaries add when the boundary conditions are the same at both of them. For 
Xi as in (|3.6|) . the asymmetry vanishes. However, the contribution of the spectral 
asymmetry will be shown to be crucial when evaluating the finite-temperature 
fermion number, which will be done in section |SJ 



3.2 6 = 7i 

Also in this case, there are no zero modes. 

At variance from the case 8 = 0, no solution exists for k = 0, and only the 
symmetric part of the spectrum appears. The eigenvalues are given by 



w n>! = ±y ^(n+ -)-J +A; 2 , with n = 0,...,oo. (3.11) 

The absence of a nonsymmetric spectrum is also easy to understand from our 
result in [H] , where it was shown that the sign of the boundary contribution to 
the asymmetry changes under a change of the intermediate sign in the projector 
defining the boundary conditions. So, in this case, the contributions from both 
boundaries cancel each other. 



4 Free energy 

With the eigenvalues of the Euclidean Dirac operator at hand, we can now 
obtain the partition function, which is 



\ogZ = log det{ipBc) 



(4.1) 



When the determinant is defined through a zeta function regularization |12| . 
one has 



\ gZ = - — 
as 



IPBC , 

a 



(4.2) 



Here, ((s, t ^ BO ) is the zeta function of the operator t ^ BO [TB| and, as usual, 
a is a parameter with dimensions of mass, introduced to render the zeta function 
dimensionless. 



4.1 6 = 

We will first discuss in detail the case where the boundary conditions are de- 
termined by 9 = 0. We will then have two types of contributions. The first 
one comes from the symmetric part of the spectrum, equation (|3.9() : the second, 
from the "nonsymmetric" part, equation i|3.10|) . These two contributions are 
given by 



6 



1 ds 



Ci(*) : 



where 



Ci( s ) = (i + (-i)- s ) £ 



i 

-00 



/ 717T \ ; 

\aLJ 



(21 



(4.3) 



(4.4) 



and 



where 



2 (is 



C 2 (s). 



(4.5) 



s=0 



(2Z + 1) — 



&(-) = E 

= (1 + (-!)-) ^ 



(21 + 1) 



a/3 



(4.6) 



Note that, in equations l|4.4|l and l|4.6[) . (— l)~ s is undetermined (except at 
s = 0), since it can be taken to be e ±J7rs , depending on the election of the cut 
when defining the complex power. As a consequence, the determinant could get 
an undetermined phase |17II18| . We will come back to this point later on, when 
performing the s-derivatives of both zeta functions. 

In order to do so, we must first perform an analytic extension of both ex- 
pressions, which are convergent for 3?(s) big enough. We first do it for £i (a 
general discussion of the procedure to be employed can be found, for instance, 
in [inj). 

We write equation (|4.4|) as a Mellin transform 



&(*) 



(l + (-l)- s ) 



(4.7) 



n = 1 

I = — oo 



Now, using the definition of the Jacobi theta function 



Q 3 (z,x) = ^2 e ~ 



(4.8) 



(i can be rewritten as 



CiOO 



(i + (-ir s ) 

-2t At 



E 



03 



(a/3) 2 ' (a/?) s 



(4.9) 
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To proceed, we will use the inversion formula for the Jacobi function 



(4.10) 



together with the definition (|4.8I) . thus getting 
(l + (-l)- s )/3a " 



Ci(«) 



2(V^) s r(§) 



E 



i-i e --(^) 2 



AttI , 



1=1 



(4.11) 



Now, the integrals can be performed to get 



&00 



(1 + (-!)-)/? 
2a-(V^) s r(|) 

'8L 



l\ 7T 2 



E 



_ s Cfl(s-l)+ 
, nlir8 

K a 



(4.12) 



This completes the analytic extension of £i. Note in particular that, due 
to the behavior of the Bessel function K v (z) for large \z\ the series in the 
second term between square brackets converges for all s. 

We now extend £2- This is quite a simple task. In fact, from the definition 
of the Hurwitz zeta function one can rewrite equation l|4.6[) as follows 



8a 



(4.13) 



We can now go to the evaluation of both contributions (Ai and A 2 ) to the 
logarithm of the determinant. Before actually doing so, note that £i and C2 
vanish for s = 0. This renders the ambiguity in defining (— l)~ s irrelevant, 
even when taking the s-derivative. For the same reason, all dependence on the 
unphysical parameter a disappears from such derivative. So, it is easy to see, 
from equations l|4.2|l . 1)4.3(1 . (|4.5|) . H4.12|) and (|4.13|) . that 



logZ 



8 


7T 


4 


2 






/9 


7T 


4 


2 


6i 





(-1)' _^ 
__ e , 



n,Z=l 
00 

log (1 + e~ 

n=l 

From this result, the free energy can be obtained 

Tog(l + e 



log 2 
hlog2. 



(4.14) 



1, 7T 2 

-logZ = 

8 B 12L 8^ 

71=1 



(4.15) 



8 



It is easy to see that, in the limit /3 — > oo, the right result for the Casimir 
energy at zero temperature (equation (|2.9(l ) is obtained. 

The high temperature (f3 — > 0) behavior of the free energy can be obtained 
by using the Euler-Maclaurin expansion, thus getting 



F = 



7TL 



(4.16) 



Note the contribution from (2 cancels, in this limit, a term linear in the 
temperature appearing in the Euler-Maclaurin expansion. 

4.2 6 = it 

In this case, there is only one contribution to the zeta function and, consequently, 
to the partition function, i.e. 



C(s) = (i + (-ir s ) ]T 

n = 

/ = — 00 



aL 



(2Z + 1) — 



(4.17) 



The steps leading to its analytic extension, and to the evaluation of the 
partition function are essentially the same as in the previous subsection. As a 
result, one obtains 



as) 



(l + (-l)-*)/3 
2a-(0F)*r(§) 



S — 1 \ 7T 2 



L 1 



1 



— Cff [s - 1, - 



and 



n = 
I = 1 



i°^=^+2f> g( i + e-^) 



(4.18) 



(4.19) 



Thus, 



24L 



^log(l + e 

n=0 



(4.20) 



In the /3 — * 00 limit, one reobtains the Casimir energy in section|2] The high 
temperature limit coincides with (|4.16|) . 
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5 Fermion number 



To evaluate the finite temperature fermion number, we will introduce a small 
chemical potential (the meaning of "small" will be specified, for each value of 9, 
in the corresponding subsection) . Then, the finite temperature fermion number 
will be calculated as 

where Z is the grand-canonical partition function. In the language of ther- 
mostatistics, this is the mean particle number of the Fermi-Dirac gas in the 
grand-canonical ensemble. The answer to the question posed in reference [H] 
will be found from the value of this object at /i = which includes, for each 
value of 9, both the zero-temperature fermion number (given in equation l|2.10|l 
and the paragraph following it) and its temperature-dependent part. 

We will follow [20] in introducing the chemical potential as an imaginary 
Aq = ijj, in Euclidean space (or, equivalently, a real A$ in Minkowski space- 
time). The Dirac eigenvalue equation then becomes 

i l0 ^ = ^ , (5.2) 

while again satisfies antiperiodic boundary conditions in the "time" direction, 
and the ones given by (|3.3|l in the "space" direction. 

The chemical potential can be eliminated from the differential equation 
through the transformation 

$ = er^oy . (5.3) 

Thus, one gets for \f the same differential equation as in the free case. 
Moreover, Vl" satisfies the same "spatial" boundary conditions but, in the "time" 
direction, one has 

#'(/3,iei) = -e^'(0,xi). (5.4) 

So, the effect of the chemical potential is to replace the "temporal" eigen- 
values in equation with 

7T 

A = A — ifi = (21 + 1)— — ifj,, I = — oo, ...,oo . (5-5) 

P 

5.1 = 

As before, we must consider two contributions to log Z 

d 



1 ds 



CiOO, (5.6) 

s=0 
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where 

oo 

Cl ( s ) = (! + (-!)-) 

n = 1 
I = — oc 



/ 717T \ z / . , „ . 7T U 

\aLJ \ ap a 



(5.7) 



and 



where 



2 ds 



Ca(s): 



s=0 



(5.8) 



&(-) = E 



Z= — oo 



ap a 



(5.9) 



In order to perform the analytic extension of d, we will proceed much the 
same way as in the /i = case. However, in order to properly write 15.71) in 
terms of its Mellin transform, and freely interchange integrals and double series, 
we will take /i to satisfy \fi\ < ~. Thus, we can write 



(1 + (-!)-) 

r(§) 



n = 1 

I = — oo 



e -4(s) L +((2'+i)^-^) L j (5 10) 



This can also be written as 



CiW = 



i Jo 



(l + (-l)- s ) 

-It ( 1 in \ At 



e 



a/3 \a(3 an J ' (a/3) 



(5.11) 



From here on, the same steps as in the previous section can be followed, to 
obtain 



Ci(s) 



S — 1 \ 7T 



f3L 
2 



E (-1)' 

n,l=l 



cosh (pdl)Ks-i 



nln[3 
L 



(5-12) 



Note the double sum in the last term between square brackets is convergent 
in the range of /i considered. So, again, d(0)=0 
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The analytic extension of £2 requires, in this case, a more careful treatment 
than in the case fi—0. In fact, 



C2O) 



E 

/— — 00 

'271-' 

a/3 

2tt 

a/3 



ap a 



E 

1=0 



-1 — S OO r 



E 

i=0 



-(H — ) — * — 



Ch[s 



"2V 



E 

1=0 



(5.13) 



Now, in order to write the second term as a Hurwitz zeta, we must relate 
the eigenvalues with negative real part to those with positive one without, in 
so doing, going through zeros in the argument of the power. Otherwise stated, 
we must select a cut in the complex u> plane |17|. This requirement determines 
a definite value of (— l)~ s , i.e., (— l)~ s = e ^ sl 9n(p-) s _ Taking this into account, 
we finally have 



pa 



1 iufi 



2tt 



(5.14) 



From (|5.12() and l|5.14|l both contributions to log Z can be obtained. They 
are given by 



a 00 

-^L + ^log(l + e-^ +2cosh(^)e- 1 ^) (5.15) 



and 



2tt 



' n(3 

log 2 + log cosh I — 



2tt 



2 2tt / 
(5.16) 



Both expressions can be seen to reduce to the corresponding ones in the 
previous section when [i = 0. 

Putting both pieces together, we finally have 

logZ = --gL +^log(l + e -^ +2cosh( A/ /3)e- :i ^) 



12L 

n— 1 

log 2 + log cosh 



M/3 
2 



(5.17) 



From this result, we can evaluate the free energy 



F 



7T 1 

12L ~ 



^log(l + e" 2 ^ + 2cosh(^/3)e~^) 



log 2 + log cosh 



M/3 
2 



(5.18) 
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It is easy to see that its j3 — ► oo limit is , independently of the value of 
This is consistent with the fact that the chemical potential has been introduced 
as a purely imaginary Aq gauge potential in Euclidean space. This corresponds 
to a real Aq potential in Minkowski space-time. Now, as is well known, this last 
can be eliminated at zero temperature through a gauge transformation. 
The high-temperature limit of the free energy is — &i — + % ■ 
According to l|5.1|) . the finite-temperature fermion number is given by 

1 + e - I r i +^ l + e-^-^ 
(^j-\sign^)Y (5.19) 

In particular, for \x = 0, it is clearly undefined. Note that this discontinu- 
ous character of the fermion number comes from zero "spatial" eigenvalues or, 
equivalently, from the spectral asymmetry. More precisely, it originates from 
the phase of the determinant. We will comment on the interpretation of this 
result in the last section of the paper. 



N 



E 



■ tanl 



5.2 e = 7r 

Since, for this value of 9, only the symmetric part of the spectrum contributes, 
we must consider the zeta function 



c(.s) = (i + (-ir s ) J2 

n = 

I = — oo 



aL 



-;-((2Z + l)^--^ 



G5.20) 



whose analytic extension can be obtained with the same method as before, and 
is given (for < ^-) by 



C(*) = 



2a-(0F) s r(f) 



S — 1 \ 7T 2 



— CH [8-1, 



1 



E (-i)' 

n = 
I = 1 



cosh (p0l)Ks-i 



(n + |) Ztt/3 X 



L 



(.5.21) 



The grand-canonical partition function, obtained by evaluating its s-derivative 
at s = with reversed sign, gives us 



logZ = 



/?7T 

24L 



X> g (i 



e 1 + 2cosh(^/3)e 1 ) (5.22) 



n=0 
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In this case one has, for the free energy, 
F =-24L-^^ l0g ( 1 + e 



2 cosh (fj,@)e 



(5.23) 



Its low-temperature limit is — 2iL wrnc hj again, is //-independent. The high- 
temperature limit is —fm — ^T- 
The fermion number is given by 



N 



E 

n=0 

which vanishes for fi = 



+M/3 



W3 



1 + e l l + e~^^ m/3 



(5.24) 



6 Discussion of the results 

The result for the finite-temperature fermion number in the case 8 = it (equation 
H5.24f> has a clear interpretation: no fermion number is created when rising the 
temperature, while keeping \x — 0. Moreover, in this case, the low temperature 
limit of the fermion number also vanishes for ji ^ 0. 

Now, in the case — 0, the result in l|5.19|l requires a careful analysis. From 
the point of view of field theory, the right answer must be taken as one of the two 
possible limits. In fact, at zero temperature, a zero mode of the Hamiltonian 
appears in Minkowski space for [i = and, as already discussed (see section [5} 
there will be two nonequivalent vacuum states, with N = ±-|. Once the theory 
is quantized around one of these vacuum states, no extra fermion number (£7(1) 
charge) will arise at one loop, when compactifying the "temporal" coordinate 
with an antiperiodic twist. This is, for 8 = 0, the answer to the question posed 
in [Hj. However, the low-temperature limit of the fermion number can be seen 
to vanish if \x is kept different from zero, due to thermal averaging over the 
degenerate ground states. 

It is interesting to note that our result for the partition function in this 
(Ramond) case doesn't coincide with the analytic extension of the one given, 
for instance, in Chapter 10 of reference |2J (see also JOI); where the contribution 
from the phase of the determinant doesn't appear. Were one to disregard this 
term, the low-temperature limit of the fermion number would vanish for \i = 
(instead of picking one of the two possible vacuum values), while it would be 
^sign(n) for /i ^= 0. 
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